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THE  CURRENT  ALGEBRA  OF  GLOBAL  MHD  STABILITY 

by 

Daniel  R.  Wells 

Department  of  Physics 
University  of  Miami 
Coral  Gables,  Florida 


ABSTRACT 

A  theory  of  nonlinear  global  magnetohydrodynamic  stability  is 
described.  The  formalism  is  an  entirely  new  approach  to  the  problem. 
The  concepts  of  space-time  and  generalized  gauge  symmetries  of  the 
flow  fields  are  invoked  to  find  constants  of  the  motion.  The  con¬ 
stants  correspond  to  charge  operators  in  a  theory  of  the  current 
algebra  of  the  fields.  The  charges,  in  turn,  are  defined  by  integrals 
that  are  determined  by  the  symmetries  of  the  fields.  The  strengths 
of  the  individual  components  of  the  currents  determine  the  amount 
of  symmetry  breaking  in  each  physical  situation.  The  constants  of 
the  motion  correspondi ng  to  the  charge  operator  are  used  in  conjunction 
with  the  principle  of  least  constraint  to  generate  the  Euler-Lagranye 
equations  correspondi ng  to  stable  plasma  motion.  For  every  symmetry 
there  is  a  correspondi ng  conserved  integral  or  charge  (Noether's 
Theorem).  The  principle  of  least  constraint  states  that  if  the  total 
enerav  of  the  flow  field  of  a  bounded  plasma  cell  is  varied,  subject 


to  a  set  of  constraint  integrals,  then  the  fewer  the  number  of  constraint 
integrals  applied,  the  more  stable  the  resulting  flow.  The  constraint 

integrals  which  generate  a  linear  (superposable)  field  yield  a  set  of 

equations  which  describe  force-free  col  linear  flow.  If  linearity  is 

sacrificed,  then  fewer  constraints  can  be  used  and  many  other  types 

of  flow  structures  are  possible. 

The  symmetries  and  corresponding  flow  structures  are  classified  by 
the  Lie  algebra  of  the  currents  and  charges.  The  formalism  leads  to 
uniqueness  theorems  necessary  to  calculate  the  type  of  structure 
present  for  a  given  set  of  boundary  conditions.  It  is  demonstrated 
that  a  multiple  integral  variation  problem  can  be  related  to  the 
principle  of  least  constraint.  The  fundamental  variational  formula 
for  the  appropriate  tensor  fields  is  developed.  The  group  generators 
for  a  particular  space-time  and  gauge  transformation  are  then  used 
tc  demonstrate  that  the  variational  approach,  of  Woltjer  and  Wentzel 
is  a  special  case  of  the  problem  of  Lagrange  wit.,  expanded  integral 
constrai nts . 

Proof  that  the  theory  actually  describes  the  lowest  lying  and 
nost  stable  energy  states  of  the  flow  structures  is  given  in  a  succeeding 
paper  -in  the  form,  of  experimental  data  which  is  a  measure  of  the  actual 
magnetic  fields  trapped  in  the  linear  structures.  Agreement  with  theory 
is  excellent. 
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Section  I  -  Introduction 

In  two  previous  papers  (Wells  and  Norwood,  1969,  and  Wells,  1971)  hereafter 
designated  I  and  II  respectively,  a  new  method  of  calculating  the  structure  of 
naturally  occurring  stable  plasma  cells  was  outlined.  The  method  invokes  a 
variat'  lal  principle  in  which  the  total  energy  of  a  closed  (bounded)  plasma 
configuration  (cell)  is  varied  subject  to  a  see  of  constraint  integrals  on 
the  flow.  The  resulting  Euler-Lagrange  equations  describe  the  magnetic  and 
flow  fields  in  the  plasma  "bunch"  or  plasmoid.  These  differential  equations 
can  be  solved  subject  to  appropriate  boundary  conditions  to  give  a  quantitative 
description  of  the  structure  of  globally  stable  plasmoids;  i.e.,  the  stability 
is  calculated  for  the  whole  plasmoid  (plasma  cell)  in  its  actual  geometric 
configuration  and  ambient  plasma  surroundings. 

The  low-lying  cell  energy  states  are  analogous  to  the  low-lying  stable 
states  of  an  atom  in  which  the  lowest  lying  and  most  stable  atomic  configura¬ 
tions  correspond  to  minimum  energy  and  the  higher  energy  states  correspond  to 
structures  that  are  more  easily  perturbed  and  broken  down  when  subjected  to 
external  disturbances.  The  type  of  stable  cell -state  changes  as  one  changes 
the  type  and  number  of  constraint  integrals  used  in  the  variational  calculation. 

The  total  energy  is  used  in  making  the  calculations  because  nearly  all 
plasmoids  (cells)  are  interacting  and  exchanging  energy  with  their  surroundings. 
This  means  that  the  structures  are  non-conservative  systems,  and  the  usual 
"’ethodj  or  suability  calculation  utilizing  effective  potentials  have  no  meaning. 
The  method  is  interesting  because  the  calculations  include  a  consideration  of 
all  nonlinear  states,  they  are  global  and  not  local,  they  make  no  assumption 
about  the  strength  of  coupling  between  various  modes  in  the  plasma  and  include 
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the  dynamic  nonlinear  terms  in  the  equations  of  motion  or  equivalently  in  the 
equations  describing  conservation  of  formal  currents  and  charges  in  the  plasma. 

The  principle  of  least  constraint  (discussed  in  II)  is  invoked  in  order 
to  find  the  various  lowest  iying  cell  energy  states.  The  fewer  the  number  of 
constraints  applied  to  the  system  in  performing  the  variational  calculation, 
the  more  stable  the  corresponding  plasmoid.  It  was  also  shown  in  II  that  the 
lowest  lying  states  correspond,  under  certain  conditions,  to  superposable  flows 
and  fields.  All  other  states  are  non-superposable  and  nonlinear.  Two  of  these 
nonlinear  states  will  interact  to  form  new  states,  rather  than  simply  superpose 
to  form  a  composite  structure. 

We  will  rederive  the  entire  theory  utilizing  an  approach  that  is  more 
rigorous  and  physically  more  transparent  than  the  methods  utilized  in  I  and  II. 
The  methods  and  formalisms  of  field  theory  will  be  used  to  derive  an  operator 
formalism  that  will  enable  us  to  extend  the  concept  of  the  stable  plasma  states 
to  a  method  of  calculating  the  growth  of  motions  and  currents  from  an  arbitrary 
initial  plasma  state. 

There  has  been  much  discussion  about  the  validity  of  this  approach  for 
calculating  the  stability  of  plasma  confinement  schemes.  A  second  paper 
follows  this  one  in  which  some  of  the  predictions  of  the  theory  are  verified 
in  detail  in  the  laboratory. 

In  this  second  paper  (Nolting,  Jindra  and  Wells,  1972),  we  will  describe 
a  plasma  confinement  scheme  (TRISOPS)  which  utilizes  all  of  these  results  to 
guide  the  production  of  high  temperature,  high  density  plasma  with  a  long 
particle  confinement  time.  Experiments  will  be  described  that  compare  the 
structure  of  the  lowest  lying  plasma  c>-ll  states  ,1'th  the  theoretically 
predicted  structure. 
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Section  II  outlines  a  derivation  of  the  formal  charges  and  currents,  dis¬ 
cusses  the  principle  of  least  constraint  and  illustrates  why  the  lowest  lying 
linear  cell  states  have  minimum  "-free  energy"  available  to  drive  "instabilities". 
Section  III  develops  an  operator  formalism  that  can  be  used  in  a  general 
attack  on  the  problem  of  the  growth  and  dynamics  of  interacting  plasmoids  (cells). 
Section  IV  applies  the  concept  of  broken  symmetries  to  plasmoid  interactions. 
Section  II  formal  Cjv^es  and  Currents 

In  I  and  II  the  lowest  lying  stable  plasma  cell  states  are  derived.  Ihe 
constants  of  the  motion  corresponding  to  space-time  and  gauge  symmetries  are 
derived  using  the  Clebsch  potentials  and  the  concept  of  a  "generalized  gauge 
transformation".  In  this  section  we  consider  all  of  the  constants  of  the  motion 
corresponding  to  space-time  and  gauge  symmetries.  We  discuss  the  expansions 
uf  Liie  corresponding  functionals  and  utilize  them  in  a  fundamental  variational 
formula  to  show  thou  minimum  constraint  corresponds  to  minimum  total  energy. 

We  then  show  that  minimum  total  energy  corresponds  to  minimum  "free  energy" 
available  to  drive  instabilities  in  the  case  of  the  lowest  lying  linear  states. 

In  II  we  discussed  the  fundamental  importance  of  Noether's  theorem  in 
relating  the  symmetry  properties  of  the  Lagrange  density  of  the  MHD  fields  to 
the  constants  of  the  motion  for  closed  plasma  structures.  It  is  necessary  to 
give  a  precise  formulation  of  the  basic  problem  of  the  calculus  of  variations 
if  further  progress  is  to  be  made  in  understanding  the  global  stability  problem 
for  bounded  plasma  cells.  We  merely  outline  the  necessary  theorems.  A  rigorous 
treatment  of  the  problem  is  given  by  Runri  (Rund,  1966). 

Given  n  real  variables  x1 ,  together  with  m  independent  real  variables  t‘, 
(Latin  and  Greek  indices  run  from  1  to  n  and  1  to  m  respectively)  consider  the 
space  Kn+rr,  or  n-n,i  dimensions  of  x1,  ta.  Denote  3  subspacc  Cm  in  R  (. 


5 


A  set  of  n  equations  of  the  type 

i  i  / . ju 
x  =  x  (t  ) 


(i) 


defies  this  subspace.  It  is  assumed  that  we  can  form  the  derivatives 
i 

— —  or  C  .  We  denote  these  derivatives  as 

.  i  m 

J  z 

p  (»‘> .  ■ 

■l  at 

Assume  that  G.  denotes  a  fixed  simply-connected  domain  in  the  m-dimensional 

space  of  the  t',  bounded  by  a  hypersurface  jG^.  Each  point  of  this  surface 

corresponds  to  a  set  of  values  of  the  t'.  Consider  a  second  set  of 

equations  of  the  type  x1  =  x^t')  representing  another  subspace  C  of  R  ,  . 

r  m  n+m 

We  require  that  this  second  subspace  coincide  with  C  for  these  values  of  t' 

m 

which  define  the  boundary  of  G^.  Then 


x 1  ( t  ‘ )  =  x i ( t  ‘ )  =  f 1 ( t  ‘ )  for  t '  «  3 , 

where  the  functions  f^t1)  are  separately  prescribed. 

Consider  a  suitably  differentiable  function  L=L(t  *,xJ  ,xJ4)  with  m+n+mn 
arguments.  This  function  is  defined  as  a  function  of  t‘  over  each  subspace  of 
the  type  of  Eq.  ( 1 ) . 

We  now  form  the  integral 


I (Cm)  =  f  L(t',xj,ij  )d(t), 
Gt 


(2) 


where 


def 

d ( t )  =  dt1 


dt 
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The  value  of  Eq.  (2)  depends  on  C  ,  i.e.,  the  choice  of  the  functions  defined 
by  Eq.  (1)  together  with  their  derivatives.  The  fundamental  problem  is  to  find 

i  ex 

the  necessary  and  sufficient  conditions  the  x  (t  )  must  satisfy  in  order  to 
yield  an  extreme  value  of  the  integral  (2). 

In  order  to  proceed  with  the  proof  another  set  of  n  equations  of  the  type 


X1  =  x’(t‘)  u) 


is  considered.  These  represent  a  1-parameter  family  of  m-dimensional  subspaces 
Cm(n)  of  R  Consider  two  neighboring  subspaces  Cm(uQ)  and  Cm(i  )»vtere  (ih^)  is 
considered  to  be  small,  and  quantities  of  order  ju-u^J2  are  neglected.  Let  P  and 
P'  be  points  on  Cm(u  )  and  Cm(u)  respectively,  corresponding  to  the  same  ta-values 

The  components  of  the  displacement  PP '  in  Rn+m  are  given  by  (0,....0;  -x ' . . 

m,  , 

x  ),  where 


lx11 
;u-u  )  —i 
1  o'  hu  i 


(3) 


u=u. 


It  can  then  easily  be  shown  (Rund,  1965,  p.  213)  that 


c1  =  (x1) 


t 1  + 


(4) 


u=u. 


and 


*• i  _  ,  *  i , 

x  =  -  ( •  x  ). 

•  t 1 

If  one  defines  the  first,  variation  of  I  (as  given  by  (2))  as 


I  .  I(u)  -  I(u0>  -  (u-u0)  (iji)  . 

wo 


then  it  can  be  shown  that 


=  V  +  A-  sV }  d(t)  +  i  —  (  L  't')d(t) 

■x-  -,k  a  'G.dtu 

i  z 
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This  can  be  written  as 


•  I 


j  L 


Ml 


Gtlnx 

using  Eo  (4). 


dt' 


hx1 

i  o i 


0  x1d(t)  + 


■  df 

bt  ?- 


l .t  4  L,  y 'c) 


•  i 

■X 


The  first  variation  takes  the  final  form, 


I  - 


•  L 


!  L  ) 


dt1 


•x 


i 

Vd(t)  +  !  —  S  ( L 


dt" 


.  X  ;  L 


Jt  ! 


+  ~r  ^x1  id(t)  .  (6) 

■x 

1  j 

Equations  (5)  and  (6)  are  the  fundamental  variational  formulas  for  multiple 
integrals.  The  integrand  of  the  second  integral  is  formally  u  divergence. 

It  is  now  easily  shown  from  Eq.  (G)  that  an  oyTre'".; ..  ioi  the  integral  defined  by 
Eq.  (2)  is  uniquely  determined  by  the  solutions  of 


where 

'<  =  1  ,  . .  n . 

This  set  of  n  equations  reduces  to  the  Eul er-Lagrange  equations  for  single 
integrals  when  m-1 .  They  are  the  Eul er-Lagrange  equations  for  multiple  integrals. 
Written  out  in  full  they  take  the  form 


M 
•  k  •  i 


v  L 


.  <  i .  •  k 

t  <t  'X  ^x  it 


.  *  •  t 

t  X 


I 


(7) 
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It  can  be  shown  (Rund,  1966)  how  certain  invariance  properties  of  multiple 
integrals  imply  that  various  quantities  are  constant  along  an  extremal.  These 
results  can  be  formulated  in  terms  of  Noether's  Theorem.  The  so-called  con¬ 
servation  laws  then  easily  follow.  Application  of  Noether's  theorem  to  MHD 
global  stability  is  discussed  by  Wells  in  II  (Wells,  1970).  Here  we  utilize  the 
formalism  of  the  fundamental  symmetry  transformations  to  directly  tie  together 
Noether's  theorem  and  the  principle  of  least  constraint. 

Consider  an  r-parameter  Lie  group  operating  on  the  variables  (tV1).  A 
typical  element  of  the  group  is  the  transformation 

t  •  =  r  (t',xj,<s),  x1  =  x1  (t' ,xj,u$)  .  (8) 

where  *  (s,t--l,  ....,r)  represent  the  r  parameters  of  the  group.  We  assume 
that  the  identity  transformation  of  the  group  is  given  for  0.  The 
infinitesimal  transformations  corresponding  to  the  finite  transformations 
given  by  Eq.  (8)  are 


where 


1  =  -(s)'s’  'x1  =  r’(sfs 


■t'  C  t  ’  -xj,a.  ) 


(summation  over  s), 

hx(t '  ,x^  ,  >t) 


(9) 


Ms) 


1 


'Vo  ‘t=0 

Corresponding  to  the  infinitesimal  increments  given  by  Eq.  (9)  we  have, 

f rom  Eq .  (4)  , 


i  i  i  .. 

x  -  x  *  t  =  - 


,i)  xs 


(summation  over  s '  , 


where 


(s)  (s)  (s) 


-x  7  x  (summation  over  s). 
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The  variations  (9)  induce  a  variation  5l  of  the  fundamental  integral  (2) 


which  we  evaluate  according  to  (5).  This  gives 

:t)  ♦  f  i  f  E 


I  =  ! 


r 

i 

aL 

d 

faL 

.  i 

dX 

’  dta 

ai1 

aj 

w«.d(t) 


’(s)V 


!Gt  *** 


,  .  3L  i 

Hs)  >•)  "(s) 

0 


,d(t).  (10) 


(summation  over  s) 

We  now  require  that  this  fundamental  integral  be  invariant  under  (9)  up  to 
an  "independent  integral".  This  means  that  L  transforms  under  (9)  according  to 


L(tfc,xJ,xj)d(t)  =  j  fL(te,xJ,x^)  +  *}d(t), 
Gt  '  Gt 


where  ( tc,xJ ,xJ^ ,a )  is  the  integrand  of  an  independent  integral.  We  must 


now  examine  carefully  the  meaning  of  i  because  much  of  what  follows  depends 
on  this  concept. 

A  thorougn  discussion  of  the  concept  of  an  independent  integral  is  given 
by  Rund  (Rund,  1966).  For  our  purposes  it  is  merely  necessary  to  state  that 
the  integral  of  a  divergence  depends  solely  on  the  values  of  its  argument 
functions  on  the  boundary  jG^  of  the  domain  G^.  over  which  we  perform  the 
integration.  Divergences  are  not  the  only  integrands  which  have  this  property. 
There  are  a  large  class  of  integrands  giving  rise  to  this  type  of  independence. 
These  are  referred  to  as  "independent  integrals".  Integrands  which  are 
divergences  provide  a  satisfactory  and  relatively  simple  theory  which  is 
adequate  for  our  purposes.  This  approach  was  extensively  developed  by  Weyl 
(Weyl  ,  1935).  There  is  another  special  sub-class,  furnished  by  certain 
determinants,  which  yield  a  theory  for  multiple  integral  problems.  This  work 
is  commonly  associated  with  Caratheodory  (Caratheodory,  1935).  There  is  a  more 
general  theory  discussed  and  developed  by  Rund  (Rund,  1966,  p.  250)  that 
includes  the  theories  of  Weyl  and  Caratheodory  as  special  cases.  This  theory 


is  interesting  because  it  classifies  the  independent  integrals  by  expressing 
their  integrands  as  homogeneous  polynomials.  This  is  extremely  useful  in 
treating  MHD  stability  problems  for  the  case  of  broken  symmetries,  since,  as  it 
will  be  shown  below,  it  is  necessary  to  expand  the  integrands  of  the 
constraint  integrals  in  order  to  treat  these  problems. 

In  order  to  greatly  simplify  the  discussion  of  the  general  treatment  of 
symmetry  transformations  of  the  type  given  in  (9),  we  assume  now  that  the  func¬ 
tions  :  are  divergences,  i.e.,  there  exist  m  functions  (t‘  ,xJ,xJ)  such  that 


where 


.jo  .  iiki  + 


s  +  -  -  - 


if  '•  depends  on  t*  only,  — —  and  -1-1—  have  the  same  meaning. 


We  also  assume  that  t  is  linear  in  the  parameters  u 


:(sj  (t\xJ,xJJJs. 


(summation  over  s) 


Then  we  obtain 


'(sTs 


d(t)  • 


Comparison  with  (10)  yields 


Gt  dt 


_lJ  1  -  f  i_ 

.xi!--(s)-«d(t)  -  jGt  dte 


-L  i 


’(s)  ,-i  (s)|  "  ’(s)  s 


*cd(t) 


(summation  over  s) 


i’i  ”'(s) 


U(s)  +  n(s)  *  *(s)  •  ls=1’*--r)- 


(summation  over  s) 


In  II  we  have  discussed  the  concept  of  generalized  gauge  symmetries  as 
developed  by  Calkin.  We  have  related  the  corresponding  constants  of  the  motion 
to  the  linear  superposatlc  states  of  closed  plasma  structures.  We  will  now 
show  that  Eq.  (13)  ties  together  the  concepts  of  Noether's  theorem,  least 
constraint,  and  symmetry  breaking  (i.e.,  expansions  of  in  complete  sets 
of  functions  of  the  spatial  coordinates). 

From  Eq.  (11)  we  note  that  there  exists  a  system  of  m  properly  behaved 
functions  ;l  =  ;  ‘(t",xJ)  obeying  the  divergence  equation.  In  order  to  simplify 
the  discussion  and  make  the  formalism  directly  comparable  to  that  developed  in 
1 1 ,  we  now  introduce  the  following  change  of  variables 


t  -xJ(x1  ,  ^2  »x^ » t ) , 

where  t  is  the  time.  We  assume  that  there  are  N  dependent  functions  ;a  of  these 
independent  variables,  and  these,  in  turn,  replace  the  functions  x1(t'1)  used 
previously 

x^t1)  •  ,a(xJ'),  (a,b,  ...  =  1,  ...,N). 

A  detailed  discussion  of  Noether's  theorem  and  related  problems  in  this  notation 
is  given  by  Hill  (Hill,  1951).  The  reader  is  referred  to  this  paper  for  a 
rigorous  treatment  of  the  conservation  equ~  .ions  which  we  now  develop. 
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Considering  transformations  that  arise  continuously  from  the  identity 
transformation,  it  is  sufficient  to  consider  the  infinitesimal  transformations 


and 


,k 

x 


k 

x 


+  ,'X 


k 


a(x‘)  =  *a(x)  +  ^a(x). 


(14) 


The  finite  transformations  are  found  by  iteration.  It  is  then  easily  shown 
(Hill,  1951,  p.  ^58)  that  one  can  associate  a  differential  conservation  equation 
with  each  infinitesimal  symmetry  transformation  in  the  form 

.  k  •  •  ;  ,  ,  -x  ,,  t 

dx  i  ;  dji  dX 

’  1  f\ 

•x 


Then  we  can  write 


7t-  +  v-S  =  0  ,  (16) 

where  the  "formal  charge  density"  associated  with  our  symmet  y  transformation 
(14)  is 


it 


3 

it 


dL 


and  the  "current  density"  S_  is 


+  >, 

"t 


(17) 


S 


t  +  (L-x - —  dx-v,‘) 


■L 

dV', 


(18) 


Then  and  the  three  components  of  Sy  can  be  associated  with  the  ms  functions 
of  Eq.  (11).  Integrating  over  a  closed  volume  •.  with  surface  /  and  using 
Gauss's  theorem,  we  have 
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TjT  [  -d(x)  =  -  {  s  d);  . 

I 

Calkin  has  shown  (Calkin,  1965  ,  Wells,  1970)  mat  one  can  write  a 
Lagrange  density  for  an  MHO  fluid  in  the  form 


wht~e  E  is  the  electric  field  intensity,  i.  and  are  the  permeability  and 

permittivity  of  the  fluid,  v  is  the  velocity  of  the  center  of  mass  of  a  fluid 
element,  .  is  the  mass  density,  p  is  the  scalar  pressure,  >.»  n,  f  are  the 
Clebsch  potentials  defined  by 


V  +  l  "  1  B  X  P  =  -  VX  +  nV 

and  P  is  a  "polarization"  vector.  The  vector  P  is  in  reality  a  type  of 
"vector  potential"  that  defines  the  current  density  according  to 

►  *•  -  V  ►  V 

j  =  -Tjr  +  7x(Pxv)  +  ( 7  •  P )  v ,  =  v  •  P  , 

and  is  the  electric  charge  density. 

Then  P  is  defined  only  up  to  a  generalized  gauge  transformation  of  the  form 
P  *  P'  +  A,  where  ~  +  vx(Axv)  =  0  and  -  0  . 

d  L 

The  las^  two  eqjation  imply  that  is  frozen  into  the  fluid.  In  II  it  is 
shown  that  the  three  vector  fields  that  move  with  the  fluid  are  B,  z  and  v 
where 

def.  _  •  - 

z  -  .  x  ( v  +  .  ’  1  GxP). 
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Then  one  can  make  the  following  infinitesimal  transformations: 

(I)  =  0  The  resulting  finite  transformations  are  just  the  canonical 

transformations  of  the  fields. 

(II)  =  ‘  <*B 

•  *•  v  •  ► 

P  -v  p1  =  p  +  aB 

This  transformation  leads  to  conservation  of  the  integral 

| A-Bdt . 

►  ► 

(III)  '  =  Mr  2 


(IV) 


r  *  P'  =  P  +  SiZ 

This  leads  to  conservation  of  the  integral 

5.'  =  ‘  i  V 

This  leads  to  conservation  of  the  integral 


I B • vd i . 


A-vdi. 


For  given  symmetry  transformations  arising  continuously  from  the  identity 
we  can  write  the  conservation  law  (16).  The  corresponding  density  n  and  current 
S  are  found  from  (17)  and  (18)  respectively.  The  appropriate  conserved  integrals 
(constants  of  the  motion)  are  then  found  by  the  application  of  vector 
identities.  In  this  way  the  conserved  integrals  in  (II),  (III),  and  (IV)  above 
nave  been  found  (see  II). 

we  now  consider  one  of  these  transformations  in  more  detail  in  order  to 
relate  our  fundamental  Eq.  (13)  to  the  principle  of  least  constraint  and  the 
concept  of  symmetry  breaking. 

Following  Calkin  (Calk’n,  1961,  p.  88)  we  consider  again  the  transformation 


P  •  P '  =  P  +  :  *B  • 


(19) 
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The  lagrangian  L  transforms  into  L',  where 

L'(  .\  .  x1)  =  L  ( ; '  ,  ,  x1)  -  S.-xE-E  . 

To  show  that  Eq.  (19)  leaves  the  action  principle  invariant,  we  write  (see 
Calkin,  1961 ,  p.  88) 

dA 


tE  * B  -  oi»( Vi  •  V*A  +  — y  *VxA) 


in  the  form 


■t 


After  considerable  algebra,  we  find 


L '  ( '-  '  ,  ,  x1)  =  L ( y 1  ,  ,  x')  +  |2-  [ff-  (A-B)  +  V-(ExA  +  »B)] 


2  L3t 


5ince  the  Lagrangian  does  not  contain  derivatives  of  the  transformed  variable 
?,  we  have  from  Eq.  (15), 


--  (A-B)}  =  -  {7- [(Ex A)  +  4>B]}  • 

If  n  is  the  unit  nurmal  to  a  boundary  surface  V,  then 


(20) 


!  2”  + 


Is-  [v-(E*A  +  $B )]dx 


=  — ^  [  {'y~  (A-B ) }di  +  (E*A-n  +  }>B-n)dJ  , 

l 

where,  for  the  transformation  (II)  above,  we  see  that 


(21) 


=  f-  (A-B)  and  S  =  |^  {E*A  +  ^B). 
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We  see  that  a  plays  the  role  of  a  "formal  charge  density"  and  S  plays 
the  role  of  a  "formal  current  density."  Thus  one  can  discuss  the  current 
algebra  of  global  stability  if  one  can  relate  these  formal  charges  and  currents 
to  the  equilibrium  and  stability  of  the  plasma  (magnetofluid)  inside  the 
boundary  surface  z  of  Eq.  (21).  For  nested  surfaces  in  an  infinitely  con- 
ducting  fluid,  the  surface  integral  goes  to  zero,  since  v,  B,  j  and  A  are  all 
assumed  to  lie  ir,  the  surfaces  (see  discussion  of  boundary  conditions  in  II) 
and  /A-Bdi  is  a  constant  of  the  motion  (see  I,  page  27).  If  both  sides  of 
Eq.  (20)  are  multiplied  by  a  time  independent  set  of  functions  of  the  space 
coordinates  (see  I,  page  32)  arranged  in  a  convergent  series,  then  where 

def. 

=  {f(x  )}a  . 

Then 

tt  H7  =  ft  i  t f(/k)  (A'B)d;  (22) 


:f(xk)  (A-B)d-r 


f(xk){7-  [(14)  +  $B] td-r  =  -  ~  f(xk)[(ExA)+$B]*ndZ 


We  assume  that  A,  B,  and  v  lie  on  nested  surfaces  of  constant  o,  p,  and  y. 
Thus  v.  ,  vp,  and  7;  are  assumed  normal  to  those  surfaces.  Since 


E  =  -  v«B 


3A 

Dt  " 


■  A 

Y  must  be  coll  inear  with  ve  if  surfaces  of  constant  are  to  te  coincident 

•  ►  * 

with  surfaces  of  constant  c,  p,  A,  B  and  v  and  if  the  gauge  symmetry  is  exact 
(not  broken).  After  equilibrium  is  achieved,  9°es  t0  zero.  The  requirement 
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A  r 

that  be  normal  to  the  surfaces  of  constant  A  is  relaxed  in  the  case  of 

broken  symmetries  since  then  v*B  is  no  longer  normal  to  the  nested  surfaces 

(there  is  leakage  of  magnetic  field  and  mass  through  the  surface  of  the  plasma 

s :ructure) . 

Returning  now  to  Eq.  (13)  and  Noether's  theorem,  we  consider  a  special 
group  of  transformations  defined  by  Eq.  (9).  If  we  make  just  the  gauge  trans¬ 
formations  corresponding  to  d,"  =  daB,  we  have  from  Eqs.  (21)  and  (22) 


—r  |od:  =  -  lr  \  f(xk)[(E*A)  +  j>B].ndr  .  (23) 

-Li  C  f 

1/ 

If  f(x  )  is  a  polynomial  with  constant  term  equal  to  unity,  then  one  term 
of  the  integrand  in  Eq.  (23)  is  (A*B).  If  only  this  single  term  in  the  integrand 
is  retained,  we  will  say  that  the  gauge  symmetry  corresponding  to  6t  =  is  exact. 
The  right  hand  terms  in  (23)  are  then  zero  and  c  is  a  constant  of  the  motion. 

If  other  terms  in  f(xk)  are  included  in  the  expansion  of  the  integrand,  we  will 
say  that  the  symmetry  is  broken.  As  we  allow  more  terms  in  the  expansion  of 

k 

f(x  ),  we  put  more  constraints  on  the  system.  The  more  terms  other  than  the 
constant  term  that  are  included,  the  more  badly  broken  the  symmetry.  The  right 
hand  terms  in  Eq.  (23)  are  no  longer  zero  but  determine  the  coefficients  in  the 
expansion  in  the  following  way. 

Let  I  =  _'f(xk){A-B}d: 


tnen 


and 


~  if (x^) ( A - B ) }dr  =  -  j  f(x^){ExA-n  +  :B*n}dl. 
f(xk) 


't  ;[1  + - Ap(c  1  MA'Bldt  =  r  ( t ) ,  wherefCxf^)  is  a  polynomial  in  xk 


and  ~(t)  corresponds  to  the  initial  value  of  the  surface  integral. 
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Now 


f  f  d’  =  3t  |  -  \  rt-"dJ 

for  convective  flow  and  any  space-time  function  F  (Wells  and  Norwood  1969,  p.  27). 
Therefore 

j^_  [{i  +  — An (^)}{A-B}]dr=  ^  [{1  +  — An($)HA*B}]di 

-  [il  +  — Afi(oi)}{A-B}]vndr. 

Therefore 

'[-i  +  — An($)HA.B>]di  -  n(t)  -  |[{1  +  —  An(«)HA*B}]v-nd£, 
or 

^  [[(Mn(:))(A-B)]dr  =  r; ( t )  -  <j  [Mn($)(A*B)]v*nd£. 

This  gives 

•  |Mn(?)(A-B)d:  =  -  ^Mn((?)«:(E<A)  +  t>B  +  (A-B)v } -ndtAt  (24) 

f  -*■  --f 

We  can  solve  Eq.  (24)  for  a  |Mn(A*B)d-.  The  right  hand  side  depends  on 
initial  and  boundary  conditions. 

Knowing  a ! M  ( )A-Bdr  for  a  give  At,  we  can  calculate  a  corresponding  Aq($) 
(Wells  and  Norwood  1969,  p.  44)  or  {An(;),  Ap(c),  Aw(c),  Ai: ( 4 )}  and  can  then 

-*■  -v 

insert  them  into  Eqs.  A(30)  through  A(33)  of  Wells  and  Norwood  to  find  aB,  aj, 

-+ 

.+  A  ft 

a v,  P,  and  ..  for  a  given  At.  This  determines  ,  etc.,  for  a  given  set  of 
initial  and  boundary  conditions  and  hence  gives  the  explicit  growth-rates  toward 
stability  of  the  fields  for  given  initial  conditions,  boundary  conditions,  and 
degree  of  symmetry-breaking.  The  method  is  independent  of  the  conventional 
linearization  restrictions  and  seems  quite  tractable  for  a  numerical  study. 
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It  is  convenient  to  define  M(tf)  in  some  applications  as  a  polynomial  with 
constant  term  unity.  In  other  applications  M($>)  must  be  expanded  in  some 
convenient  complete  set  of  functions  of  the  electric  scalar  potential  (see  II, 
p.  648). 

If  one  does  not  want  to  perform  the  iteration  calculation  indicated  by 
Eq.  (24)  but.  instead,  wants  to  calculate  the  final  cell  configurations 
that  will  result  for  a  given  degree  of  symmetry  breaking,  then  a  modified 
constant  of  the  motion  must  be  found.  In  this  case 
def 

a  =  MnU)a  » 

where  M  ( c )  is  an  appropriate  complete  set  of  functions  of  the  electric  potential, 
$• 

Then  we  obtain 

r  .  -+  -*  6a  r  •* 

~  ad  i  =  4j  i  ~  M(^)(A.B)d1  =  -  M  ( 4>){V*[(ExA)  +  *B]}dT. 

at  J  at  J  c  n  c  j  n 


This  gives 


it  | 3dl  = 


vMn(<^)' 


[(E*A)  +  <j>B]di 


'Mn(<|>)[(ExA)  + 


ij)B]di 


-  j vMn ( <p ) -[(ExA)  +  i?5jdi  -  j{Mn($)[(ExA)  +  ^B] } * nd:;  . 

£ 

One  must  now  assume  that  during  the  time  that  the  cell  decays  to  its  final 

lowest  energy  state,  the  time  average  fields  (Ex/,)  and  <t-B  remain  on  the  nested 

surfaces.  Then  for  the  boundary  surface, 

M  (e)  =  {M  ($.)}  =  a  fixed  constant, 

n'  n  £ 

Thus 


3 

3t 


ad ; 


7M^( ij) )  •  [ ( ExA)  +  $B]di  - 


<Mn(' 


[(ExA)  +  4'B] )  *nd:;  (25) 


The  first  term  on  the  right  hand  side  drops  out  because  v Mn ( 4 )  is  normal  to 

-¥  -V  * 

the  nested  surfaces.  The  second  term  on  the  right  is  zero  because  (ExA)  and  iB 
1  ie  on  the  surfaces . 
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Now  consider  a  different  group  of  transformations  of  the  type  defined  by 
Eq.  (9)  which  include  not  only  the  gauge  transformation  corresponding  to,  say, 

=  ‘‘mB  but  also  a  simultaneous  transformation  of  the  time.  Then  by  Eq.  (9) 

->  y 

t  P  P  •  > 

t  =  .  and  :x  =  cT1  i!  (summation  over  v  and  ;i)  (26) 

where 

*  ~y 

_P  p  ->  "*■  '* 

t  =  t  +  .  and  x  =x  +n=P+  6aB. 


One  must  examine  the  meaning  of  this  transformation.  Consider  the  change  in  a 
function  F(x)  under  an  infinitesimal  transformation  of  the  form 


x'  =  f  ^  ( x-j  •••  xn;  a^  •••  ar),  i  =  1,  -  ,  n,  (Hammermesh,  p.  296). 


dF  =  z  dxi  =  I  l  Uu  (x)'a> 

i  =  l  xi  1  i=l  xi  t=l  1?- 


n 


where 


a  ■  (x)  = 


k=l 


-fil 

— j  ia,. 

■’ak!  k 
a=0 


Then 


r  n  r 

dE  =  _  a  {  l  a .  ,(x)  F  =  7  5a,x„F. 

.=1  1  i=l  ^  ,xi  t=l  '  * 


The  operators 


n 

a  ■  (x)  - — 

1  =  1  1  'xi 


are  the  group  generators. 

In  the  particular  case  of  the  group  described  by  Eq.  (26),  one  has 


dP  =  i 
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We  have,  therefore 


j 

leu  1 

Z6B 
V5Y 

jne  can  now  easily  anticipate  from  this  relationship,  derived  from  the 
properties  of  the  group  generators,  the  form  of  the  Euler-Lagrange  equations 
of  the  corresponding  flow. 

For  a  nonmagnetized  medium  moving  with  a  velocity  v  which  is  small  compared 
with  the  velocity  of  light 

-* 

vxB  =  ~  +  vx(Pxv)  +  v(7-P). 

Then 

p  ' 

~  =  7XB  -  [7x(Pxv)  +  V ( 7 * P ) 3 . 

The  requirement  on  the  group  generators  is  that  ^  be  parallel  to  B,  z,  and  v. 

Then 

—■  =  iB  “  { 7xB  -  [vx(Pxv)  +  v(?*P)]}, 

=  .i  =  j { 7xv  +  -  BxP}  =  [vxB  -  vx(Pxv)  +  v(v*P)], 

,p  .  •  ■*  -*  -+  -► 

~-t  =  v  =  [vxB  vx(Pxv)  +  v(v.P)3 

where  t,  B,  and  >  are  scalar  constants. 

These  equations  are  trivially  satisfied  for  a  force-free  col  linear  flow 
on  closed  nested  surfaces.  This  is  the  flow  described  by  the  Euler-Lagrange 
equations  resulting  from  a  variation  of  the  total  energy  of  a  closed  plasma 
volume  subject  to  the  constants  of  the  motion  correspond! ng  to  these  three  symmetries. 
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The  coefficients  of  the  group  generators  then  become 


—  =  1,  and  cP 

3  v  I 


lax  _  -1 

i  —  ~  i » 


We  now  write  Eq .  (13)  in  the  form 


’dt~M  'aX1 


-L  i 

TT!  ’<*> 


*  B  £ 

-  - —  i a,  \  +  \  (summation  over  (s)),  (27) 

dt6  [  (s}  (s). 


where 


(s)  -  iv 


From  the  definition  of  r'^sy  we  ^dve 


,/  -  .  i  l  ~  V  ■* 

(S)  |K  ;i  Xa  "(s) 


I 

T  "(s)  (summation  over  s). 


It  is  now  convenient  to  define  a  Hamiltonian  tensor  (Rund,  1966,  p.  240)  in 
the  form 


H  ‘  (t  ,x^,Pj) 


-  L(t',xj,^'  (tt:,xh,pp)  +  pf  :i(tE,xh,ph) 


This  implies 


>  (s)  “  H  (s) 
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From  the  trans formations  Eq.  (26)  and  the  definition  of  njsj  given  on  page  9 


above 


or 


and 


P  _  P  p  t 

i  -  =  c  -►  -  Pf  S 
(6aB)  (6aB)  (v) 


P  _  ,'f*  p  A 

n!l  -  h  -  Pt  f’v 


1  - 


r 


vJe  then  obtain 


„6  _  u6rt  oL  P 

e(s)  ■  Ht«(s)  *  TT  Hs) 
which  reduces  to  ° 


(Summation  over  s) 


8(s)  -  Ht  -  J-  ■  "5 


6 


since  L  does  not  contain  terms  in  P 


6' 


We  have  assumed  that  ~  <5v  *  B6a 

7L 


The  left-hand  side  of  Eq.  (27)  is  zero  on  any  extremal  subspace  (i.e., 
the  Euler-Lagrange  equations  are  obeyed).  Thus,  it  takes  the  form 

-  {H®  +  =  0. 

dtH  r  n 

Utilizing  the  Lagrange  density  given  above,  the  Hamiltonian  density,  H, 
is  proportional  to  the  total  energy  per  unit  volume  (Rund,  1966,  page  300)  of  the 
fluid;  i.e.,  we  have 


We  may  write 


=  - <5  where  t4=  t . 


dt 


where 


{-  f  *  4>)  ~  0  , 

def . 


t  =  energy  density  of  the  plasmoid 

-  ,  +  ;  =  f(t;  ,t-  ,t  )  =  f(x,y,z) 

!  i ^  Id i  =  constant 
»  0 

-  E  +  j4>^j  dr  =  constant  (summation  over  s) 


def. 


;’(s) 


(S) 


the  total  energy  of  the  closed  plasma  structure 
d;  =  E  +  constant  (summation  over  s) 
f($)di  =  E  +  constant. 
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From  the  discussion  above,  it  follows  that  the  smaller  the  number  of  terms  in¬ 
cluded  in  f ( ; ) ,  the  polynomial  expansion  of  the  constraint  integrand,  the  lower 
the  total  energy  of  the  closed  plasma  structure  under  consideration. 

The  arguments  presented  above  apply  specifically  to  the 
generalized  gauge  symmetries.  For  the  space-time  symmetries,  the  Lagrange 
density  is  invariant.  Thus  the  independent  integrals  t^are  identically  zero  for 
the^e  transformations.  The  constants  of  the  motion  corresponding  to  the  latter 
symmetries  are  also  independent  integrals,  however,  and  can  be  written  as 
divergences  in  the  sense  of  Eq.  (15).  If  we  perform  the  space-time  trans¬ 
formations,  Eq.  (27)  can  be  shown  (Rund,  1966,  p.  296  and  297)  to  take  the 
form 

dH 

-r  =  0. 

dt 


Any  other  independent  integral  of  the  form 


=  0 


can  be  added  to 


dH" 

— —  =  0  ' 
dt° 


Since  the  constants  of  the  motion  corresponding  to  the  space-time  symmetries 
are  of  this  form,  the  relationship  between  the  total  energy  and  the  independent 
integrals  is  of  the  same  type  for  these  symmetries  and  the  generalized  gauge 
s ymmetri es  discussed  above.  We  have  here  an  alternate  proof  of  the  principle 
of  least  constraint  (see  II,  page  663)  i_f  the  total  energy  of  the  structure 
can  be  directly  related  to  the  free  energy  available  to  drive  various  types  of 


i 

! 


i 


2  5 


instabilities.  This  relationship  must  be  established  independently  for  each 
type  of  plasma  structure  formed  for  the  various  possible  combinations  of  con¬ 
straint  integrals  (see  I,  page  45).  In  I  and  II  we  have  described  the  force- 
free,  collinear  stricture  as  a  very  low  energy  and  very  stable  structure.  This 
structure  corresponds  to  exact  gauge  symmetry  and  totally  broken  space-time 
symmetry,  i.e.  the  constants  of  the  motions  are 

f'A-vd!  , 

B  *  vd  ;  , 
and 

A-0d: . 

Linear  momentum,  angular  momentum  and  total  energy  are  not  conserved.  The 
resulting  Eul er-Lagrange  equations  are  found  to  be 

vxB  = 
and 

v  =  j  B , 

where  •  and  are  scalar  constants. 

If  we  write  Euler's  equation  for  a  steady  state  fluid  in  the  form 

0  =  • (P  +  \  -V-  )  +  jxB  -  p(cxv) , 

the  force-free  collinear  equilibrium  equation  becomes 

C  -  v(p  +  I  .  v  ').  (2g) 

It  can  be  shown  (see  Sherd  iff,  1965 ,  p.  187)  that  the  complete  thermo¬ 
dynamic  energy  equation  can  be  written 
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j  i,  OY  (U  +  y~)  -  E*j  -  div  Q  -  div  pv  +  viscous  terms  id t  =  0, 


where 


r  '  e 

t  !  i 


E-ju  is  the  rate  at  which  electrodynamic  energy  is  supplied  to  any 


instabilities  growing  in  the  structure, 
driving  i ns tabi lit.  ;, 


divQdi  is  the  heat  flux  available  for 


div  pv  at  is  the  flow  work  available  for  driving  instabil¬ 
ities  and  the  viscous  term  is  the  energy  available  for  resistive  mode  instabilities. 
From  Eq.  (28)  and  the  relation 

U  =  |  p  , 


where  U  is  the  internal  energy  of  the  fluid,  we  have 
•(§  U  +  i  .  V-  )  =  0 


(29) 


for  the  force-free  collinear  structure.  But  the  rate  of  change  of  the  energy 
available  to  drive  all  the  instabilities  (the  free  energy  for  our  plasma  model) 
is 

jjffu*  •?•'>'&<«*  (f)  ♦  ?•<•<«  *  \  <•*•>. 

For  a  steady  state  force- free  collinear  structure  , 

Tt  <f  ^  r>  ■ 


This  implies 


Dt(|  U  +•  2  )  l'+  2 


1 


=  0, 


where  we  have  used  Eq.  (29).  We  see  that  for  this  model,  the  lowest  lying 
linear  plasma  state  has  no  free  energy  available  to  drive  instabilities. 
Rewriting  Eq.  (27)  in  the  form 


d 

dt 


'(s)  "  7s) 


)  =  0 
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*or  an  extremal  subspace,  we  see  that  the  exact  form  that  the  integrands  will 
take  depends  not  only  on  the  amount  of  symmetry  breaking  present  in  a  given 
physical  situation,  but  also  on  the  particular  representation  of  the  independent 
vecto”  fields  that  is  utilized  in  writing  the  Lagrange  density  and  making  the  Lie 
group  of  transforn,  ions  represented  by  Eqs.(9).  The  number  of  terms  retained  in 
the  expansion  f(f.)  is  determined  by  how  badly  the  assumption  underlying  the 
symmetry  is  violated.  For  example,  if  there  is  some  ohmic  resistance  present 
in  the  fluid,  then  the  lines  of  magnetic  induction  will  slip  slowly  through 
the  fluid  and  the  symmetry  represented  by  the  transformation  &  A  =  will  be 
weakly  broken.  This  amount  of  line  slip  must  be  represented  by  the  terms  re¬ 
tained  in  4  f(4)  A-Bd :  .  If  there  is  no  slip,  i.e.,  an  ideal  fluid,  the  symmetry  is 
exact  and  f(:)  =  1,  etc.  If  one  uses  a  special  representation  of  the  fields, 
for  example,  the  representation  for  azimuthally  symmetric  fields  described  in 
II,  (Wells,  1970,  p.  653),  then  the  group  parameters  c<s  in  Eqs.  (9)  will  be 
different  and  we  have  an  entirely  different  representation  of  the  same  Lie 
group.  This  will,  in  turn,  mean  that  the  integrands  in  the  constraints  will 
change  and  there  will  not  be  a  one  to  one  correspondence  in  the  integrands  cor¬ 
responding  to  the  same  degree  of  symmetry  breaking.  This  is  why  we  have  written 
in  II,  page  6b5,  Eq.  (32), 

I21  -  T(_ ■ P ) ' d n  ~  jA-tfdi 

i.e,  £-3d  replaces  T(_P)'d'  but  the  mapping  is  not  one  to  one.  The  resulting 
Euler-Lagrange  equations  are  the  same,  however,  (compare  Eq.  (38)  page  656  in  II 
witn  Eg.  (14),  page  29  of  I).  This  must  be  so  since  the  two  group  representa¬ 
tions  have  the  same  group  algebra.  It  is  interesting  to  note  that  the  expansion 

utilized  by  Woltjer  for  azimuthally  symmetric  systems,  i.e.  l  a n(wP)  ,  meets 

n=0  n 

our  requirement  that  7f(-)  =  a  (^P)n  )  is  normal  to  the  nested  surfaces. 

n=0  n 
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We  have  snown  that  Eq.  (27)  can  be  written  in  the  form 


+  =  constant  (summation  over  s)  v3C 

where  is  the  total  energy  density  and  4^ y is  the  integrand  of  an  independent 

integral  represent  q  the  sum  of  the  constraints.  These  integrands  have  been 
assumed  to  be  linear  in  the  parameters,  n.  .  Indeed,  for  our  special  problem, 
we  found  that  the  integrands  i^s ^  took  the  special  form 

(A-B),  (8*v)  and  — _  (A-v) 

-  ► 

where  ;,  and  •  are  proportional  to  the  group  parameters  isB,  m?z  and 


If  Eq.  (27)  is  considered  to  represent  a  relationship  between  the  total 
energy  of  a  closed  plasma  configuration  and  the  appropriate  constraints  on  the 
flow,  then  one  can  vary  the  integrands  and  replace  the  coefficients  etc., 
with  a  new  set  of  parameters ,  the  Lagrange  multipliers.  These  will  now  be 
uroportional  to  the  original  group  parameters.  The  old  parameters  ,  and  the 
first  Lie  group  of  transformations  transform  the  arbitrary  vector  fields  x1 ,  or 
equivalently  .  \  into  a  subspace  which  represents  equilibrium  flow,  i.e.,  the 
Euler-Lagrange  equations  are  satisfied  and  the  left  hand  side  of  Eq.  (27)  is 
identically  zero.  If  we  write  the  equation  for  this  second  variation,  Eq.  (30) 
takes  the  form, 

-  -  d ■  -  !  '(s)*(s)dT  =  0  (summation  over  s)  (31) 

lor  the  ease  of  the  generalized  gauge  symnetries,  this  equation  is  in  the  form 

of  Eq .  ( 2 i )  of  II,  i.e. 


«.  ‘I.  =  0 

in  in 


i  - 1  n=o 
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For  the  case  of  space  time  symmetries ,  is  zero,  i.e.,  the  Lagrange  density 

is  invariant  with  respect  to  the  space  time  transformations  (Calkin,  1961). 

Thus,  if  constants  of  the  motion  corresponding  to  space  time  symmetries  are 
to  be  incorporated  into  the  variations,  they  must  be  added  separately  to  Eq. 

(32).  As  we  obse.  ed  above,  however,  this  does  not  affect  any  of  the  present 
arguments . 

Equation  (31)  is  a  special  case  of  "the  problem  of  Lagrange"  (Rund,  1966, 

P.  323).  In  this  problem, the  curves  which  afford  extreme  values  to  the  funda¬ 
mental  integral  are  required  to  satisfy  certain  subsidiary  conditions.  These 
may  be  in  terms  of  first  order  differential  equations  or,  in  our  case,  constraint 
integrals . 

The  Lagrange  multipliers  ^are  a  new  set  of  parameters  correspondi ng  to 
a  second  Lie  group  of  transformations  that  take  the  equilibrium  vector  fields 
of  the  first  variation  and  transform  them  into  the  special  subset  of  fields 
describing  the  equilibria  corresponding  to  the  constrained  fields.  These  may 
or  may  not  be  relatively  stable,  depending  on  the  number  of  constraints 
employed.  The  variational  calculation  of  Woltjer  and  Chandrasekhar  (See  I  and 
II)  is  simply  a  method  of  incorporating  the  Lie  group  of  transformations  on 
P  and  the  time  into  another  second  group  of  transformations  of  the  equilibrium 
fields  into  constrained  equilibrium  fields. 

It  is  interesting  to  note  that  the  reason  each  integrand  in  the  con¬ 
straints  n,:'ct  be  expanded  in  an  analytic  series  if  one  is  to  recover  all 
possible  equilibria  is  now  quite  clear.  We  have  already  determined  all  possible 
equilibria  when  we  obtained  the  fundamental  variational  formulas,  i.e.  Eq.  (6) 
and  (7).  When  we  seek  constants  of  the  motion  by  performing  the  Lie  group  of 
transformations  given  by  Eq.  (9),  we  are  projecting  out  a  subgroup  of  constrained 
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equilibria.  When  we  perform  the  second  group  of  transformations  given  by  Eqs. 
(31)  and  (32),  we  make  a  second  projection  to  a  manifold  corresponding  to  a 
further  specialization  of  constraints.  When  we  again  demand  all  solutions  of 
Eqs.  (6)  and  (7),  we  must  reverse  the  whole  projection  process  by  complete  expan¬ 
sions  of  our  constr.  nt  integrals.  This  reverse  process  is  of  practical  interest 
only  for  fields  close  to  minimum  energy  when  one  is  studying  the  rate  of  change 
of  the  stable  structure  for  weakly  broken  symmetries.  We  will  illustrate  these 
concepts  in  detail  in  Section  IV. 

We  should  note  here  that  if  one  performs  the  second  variation  to  solve  our 
problem  of  Lagrange  using  many  terms  in  an  expansion  of  one  or  more  of  the 
constants  of  the  motion  (independent  integrals)  corresponding  to  the  gauge 
symmetries,  combined  with  a  constant  of  the  motion  corresponding  to  say  one 
exact  space-time  symmetry,  the  resulting  energy  may  be  higher  than  that  resulting 
from  another  situation  in  which  all  symmetries  of  whatever  type  are  nearly  exact. 
In  this  sense,  the  energy  is  not  a  monotonical ly  decreasing  function  of  the 
number  of  classes  of  constraints  employed. 

The  details  of  these  transformations  for  the  case  of  the  exact  gauge 
symmetries  is  most  easily  developed  in  terms  of  the  "formal  charge  density" 
defined  by  Eq.  (17)  and  the  corresponding  formal  charge,  Q. 

In  a  field-theoretical  formalism,  a  formal  charge  Q  is  defined  as  the 
space  integral  of  the  zeroth  component  of  a  local  four-vector  current, 

U(xo}  =  ldV3j0(x)-  (33) 

These  quantities  appear  in  discussions  of  symmetries  and  broken  symmetries 
in  quantum  field  theory  and  are  used  in  the  "current-algebraic"  approach  to 
particle  theories. 
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In  quantum  field  theory,  relations  of  the  type  given  by  Eq.  (22)  have 
rather  troublesome  convergence  properties.  There  are  other  major  difficulties 
in  quantum  theory  which  do  not  arise  here  as  long  as  the  Lagrangian  and  the 
fields  remain  unquantized. 

In  order  to  ma-  a  clear  the  relationship  of  Eq.  (22)  to  the  formal 
charges  defined  in  Eq.  (33),  we  note  that  for  our  purposes,  jQ(x)  in  Eq.  (32) 
is  defined  by  Eq.  (22).  Then  Q,  the  formal  charge  operator,  is  given  by 
;[jQ]dV3  where  j  a.  In  the  usual  formalism,  Q  is  actually  time  independent. 

In  our  case  we  have  already  proved  that 

jod(x)  =  ~  [Q]  =  0.  (34) 

Section  III.  The  Lowest  Lying  Plasmoid  States  and  the  Charge  Operator  for 
Linear  Fields 

In  Section  II  we  derived  the  formal  charge  [Q],  i.e.,  Eq.  (34).  [Q]  can  be 

used  to  perform  a  variational  calculation  of  the  flow  equilibrium  solutions  for 
the  fluid  (see  I).  This  variational  calculation  is  equivalent  to  a  Lie  group 
of  transformations  that  take  the  vector  fields  describing  a  general  MHD  flow 
field  into  a  vector  subspace  (manifold)  which  describes  linear  (superposable) 
equilibrium  flows. 

if  one  defines  a  complex  vector  field 
♦ 

=  u  +  iB  , 

then  the  transformation  is  of  the  form; 

•  '  =  f(  -,Q). 
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is  a  vector  subspace  describing  equilibrium  solutions  that  corre¬ 
spond  to  force-free. coll  inear  flow.  This  subspace  is  defined  by  an  eigenvalue 
equation  of  the  form 


where 


A  •  7  )  B)  i  <  )  B  ' 


-»•  def. 

|B>  (ii6)B, 


<  is  a  constant  (see  I),  a  is  a  scalar  defined  below,  and 
"1  =  "  M23  +  lM14 

2  =  '  M31  +  lM24 

3  =  "  M12  +  iM34‘ 

The  M  are  defined  in  terms  of  the  Maxwell  stress  tensor. 


The  electromagnetic  force  terms  are  introduced  formally  in  the  usual  way 


=  f  J  '  , 
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where 
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J  is  the  electromagnetic  four  current. 
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The  equation  Tor  f  can  be  written 


f,v  *  if23  M23  +  if31  M31  +  ifl?  M12 

+  i f  41  +  if42  M42  +  if43  M43 

or 

[f  ]  =  i  [ifpv  M  ]. 

J  \!  C  pV 

The  matrices  M  are  the  Hermitian  and  anti-Hermitian  generators  of 

U  V 

Lorentz  transformations  (Kursunoglu,  1 962 *  page  50).  They  are  introduced 
here  because  they  will  be  used  later  in  a  group  theoretical  discussion  of  the 
formal  charge  operators. 

Maxwell's  equations  take  the  form, 
f 


Equation  (35)  is  equivalent  to  the  equations  derived  in  I  and  II  for 
describing  the  fields  for  exact  gauge  symmetries,  i.e., 
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Vxv  =  k-V 


vxB  =  <B 


v  *  v  =  0 


and  i 


v-B  =  0 


where  k  and  £  are  :onstants. 
One  notes  that 


Ai  = 


f  0  0  0  11 

0  0  i  C 
0  -i  0  0 
10  0  0 


Ao  = 


and 


We  may  now  write 
v  _ 


A4 


f  0  0  -i  01 

0  0  0  1 

i  0  0  0 
0  10  0; 

0  i  0  0 

-i  0  0  0 
0  0  3  1 

i  0  0  1  0J 

10  0  0 

0  10  0 
0  0  10 
0  0  0  1 


'l  ”1  +  aZ32  +  a333  +  A494  , 


=  ±eB 
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One  can  now  differentiate  between  the  formal  charge,  which,  in  this  case, 
explicitly  di  plays  the  generalized  gauge  symmetry  (see  II)  for  a  conserved 
current,  and  a  "differential  charge  operator."  The  formal  charge  operates 
through  a  variatical  transformation  to  define  a  Lie  group  of  transformations 
of  the  fields  related  to  a  definite  symmetry.  The  "charge  operator*  is  related 
in  this  case  to  the  differential  operator  that  defines  the  vector  subspace 
corresponding  to  the  transformation. 

We  may  wr-1  ce 


which  implies 


(36) 


In  addition  to  the  force-free  conditions,  the  linear  superposable  fields  are 
collinear  (see  II). 

This  means  that 


u  =  •  =  B, 


where 


u ' 

— -  (a  scalar), 
B 


so  that  we  may  write  the  defining  equation  for  |n>  as 


•  •  =  (i  •  ,)  B. 

For  the  special  case  of  collinear  fields. 


1 B  • . 


(37) 
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The  charge  operator  must  produce  both  of  the  transformations ,  Eqs .  (36)  and 
(,37).  We  wi  '  also  require  that  our  charge  operator  conform  to  most  of  the 
other  restrictions  usually  put  onto  similar  operators  in  quantum  field  theory. 
In  most  cases  the  '  requirements  will  not  have  the  same  physical  significance 
that  they  have  in  field  theory  since  our  theory  is  classical.  We  will, 
nevertheless,  find  the  analogy  very  useful  in  our  later  discussion  of  symmetry 
breaking  given  in  Section  IV. 

The  differential  operator  G  corresponding  to  the  formal  charge  operator 


Q 


[G]  -  M/k-  (••V)]. 


Then  for 


def . 

IB:-  ^  ( i  * 3 ) 


f  B1  ) 


G  IB-  =  |B>. 


(38) 


G  operates  in  the  vector  space  defined  by  jB'5,  since  the  formal  charge  Q 
transform:,  the  space  defined  by  all  solutions  of  the  equilibrium  equation,  i.e. 

{ ■  ,  to  the  special  subspace  corresponding  to  force-free  collinear  flow,  the 

i 

i b  .  We  have  written  B  -  as  a  three-vector  in  order  to  retain  the  M^v  (Lorentz 
group  generators).  The  would  more  properly  be  written  in  covariant  form. 

5, nee  this  would  only  make  the  formalism  more  obscure,  we  sacrifice  the  elegant  co¬ 
variant  formalism  for  the  sake  of  clarification  of  the  physical  principles  involved. 
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Section  IV  Broken  Symmetries  and  Plasmoid  Interactions  (An  Example) 


One  now  proceeds  to  apply  the  complex  operator  description  given  by  Eq.  (38) 
to  the  calculation  of  the  fields  of  two  interacting  plasmoids. 

Equation  (38)  describes  the  structure  of  a  force-free  coll  inear  plasmoid. 

Two  of  these  plasmoid  will  superpose  linearly  if  =  62  and  k.  =  for  the  two 
structures.  If  or  f  or  both  inequalities  hold,  any  interaction  of 

the  two  plasmoids  will  be  nonlinear  and  a  detailed  treatment  of  the  process 
would  involve  expansion  of  all  the  Lagrange  multipliers  and  a  machine  calculation 
of  the  resulting  interaction  (Wells  and  Norwood  1969).  This  will  be  treated  in 
detail  in  a  later  paper.  We  are  content  here  with  an  approximate  treatment  which 
we  will  apply  directly  to  a  problem  of  immediate  interest,  i.e.,  the  interaction 
of  two  plasma  vortex  structures  interacting  at  the  center  of  a  magnetic  mirror. 

The  experiment  is  described  in  detail  in  paper  II  following  and  in  (Wells 
and  Norwood  1969).  Two  plasma  vortex  structures  are  fired  at  each  other  inside 
a  magnetic  mirror  by  conical  theta-pinch  guns  placed  inside  the  mirror  system 
as  shown  in  Fig.  1  of  paper  II.  It  has  been  demonstrated  elsewhere  (Wells  and 
Norwood  1969)  that  one  structure  is  corotational  and  the  other  countrarotational . 
This  means  in  the  equation  expressing  collinearity 

v  =  bB 

both  structures  do  not  have  the  same  sign  of  6,  which  implies  that  the  interaction 

will  be  nonlinear.  We  assume  that  since  it  can  be  shown  that  ►  ^  , 

where  .  is  the  plasma  frequency  and  c  the  speed  of  light;  and  the  plasma 
parameters  are  the  same  at  each  end  of  the  machine. 
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Thus  one  sees  that  although  the  vector  fields  describing  a  force-free 


col  linear  stricture  correspond!' ng  to  unbroken  generalized  gauge  symmetries 
(described  in  II)  are  linear  fields,  and  the  free  noninteracting  structure  is 
built  up  from  those  fields,  supplementary  conditions  on  interacting  structures 
produces  a  non! i nea. _ interaction. 

The  structures  meeting  at  the  center  of  the  mirror  are  not  low  3  structures 
even  though  they  are  force- free.  (See  discussion  in  Appendix  II.)  The  inter¬ 
acting  inertial  flows  and  magnetic  fields  build  up  a  magnetic  barrier  or  vail 
(Wells  1968)  which  isolates  the  tw^  structures  for  times  close  to  zero  and  thus 
justifies  an  approach  which  might  be  termed  a  "quasi-free  structure"  approach. 

One  assumes  that  the  linear  theory  that  describes  a  free  structure  applies  for 
times  close  to  zero  but  that  the  initial  interaction  produces  a  small  breaking 
of  the  generalized  gauge  symmetries.  This  means  that  .. ,  which  is  shown  to  be  a 
constant  if  the  gt uoe  symmetries  are  not  broken  (see  I  and  II),  is  now  a 
function  of  the  spatial  coordinates. 

We  have  shown, 


Tne  numerator  and  denominator  go  to  zero  at  the  same  rate  as  .•  =  •  (i!0‘)~:/  • 

If  the  gauge  symmetries  are  broken,  then  one  can  expand  the  constraint 
integrals  in  a  complete  set  of  functions  of  the  space  coordinates,  say  M  (t) 
and  to  (:)  (see  Wells  and  Norwood,  1969,  p.  43).  This  gives 

,  „M 

■  =  f(x,y,z). 

ju- 
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We  will  assume  that  experimental  observation  indicates  that  initially  the 


symmetry  breaking  is  small.  Then  k  ~  f,  i.e.,  f  is  a  slowly  varying  function  of 
the  spacial  coordinates.  In  order  to  proceed  with  the  calculation  of  the  initial 
decay  rate  of  the  interacting  structures,  it  is  necessary  to  show  that  G  in 
Eq.  (38)  is  the  time  development  operator  for  the  "quasi-free"  structures. 

In  the  quantu  -mechanical  description  of  an  energy  eigenstate  |ur(t)>,  one 
has 

HjuE(to)-  =  E  juE(t0)>. 

where 

E  =  fix. 


We  then  may  wri te 

|ur(t)'  =  (t 


=  e-ii:(t-to)/*|UE(to)> 
which  may  be  put  into  the  form 


|uE(t)-  =  e'iH(t-to)/fliuE(t0)-.  (39) 

This  defines 

U(t,tQ)  -  e'iH(t-to)/h 

where  U( t , tQ)  is  the  unitary  time  development  operator  for  the  eigenstates  1 u£ ( tQ )  - 
In  a  force-free  collinear  structure,  both  the  velocity  and  magnetic  field 
are  described  by  differential  equations  of  the  form  (Bjorgum  and  Godal  1952). 

vxc  =  <6  v*c  =  0 


V '  C  +  «':c  =  0 

If  the  medium  has  viscous  losses, the  equilibrium  equation  for  the  structure 
can  be  written 

-  ux’.’xu  =  F  -  -  vp  -  v  ~  +  V-  v2u, 

OL  CQ 

¥  ►  > 

where  F  =  jxB  =  0  (force-free)  and  -  v  *  viscosity  coefficient  . 

P 
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From  these  equations  one  has 


-^jr  +  ( v  *  *' )  u  =  0, 

which  may  be  solved  by  the  method  of  separation  of  variables 
u  =  f  (t)v'/x,y,z) 

to  give 

--  -  V  <  ^  t  -*•  /  \ 

U  =  e  w(x,y ,z) . 

This  implies  that 

;u(t)>  =  e'v<‘t|u(t0)>,  (40) 

but  from  Eq.  (38l 

A  •  7  i  B  =  i.:  jB1' 

so  that  in  complex  operator  space 

|B(t)>  =  fev^'7^(t‘toh|B(t0)>,  (41) 

where  Equations  (37)  and  (40)  were  used.  Comparison  of  Eqs.  (39)  and 
(41)  indicates  that 

U(t,tQ)  -•  exp  (v<  G2)(t-tQ). 

The  value  of  can  be  calculated  independently  and  is  a  function  of  the 
plasma  temperature  and  density.  U(t,t0)  is  a  hermitian  time  development  operator 
that  is  generated  by  (•,■•- G:).  Thus  one  has  proved  in  this  approximation  that  G;' 
is  the  generator  of  the  time  development  operator  for  the  decaying  structures. 
Infinitesimal  unitary  transformations  may  be  written 

U  =  1  +  i- 0G, 

where  is  infinitesimal  and  0r  is  Hermitian. 
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An  operator  -i  changes  under  an  infinitesimal  unitary  transformation  in  the 
following  way: 

r  =  =  (1  +  if0gji:(l  -  icO^)  =  s:  +  i t  [0o »  _]  (to  first  order  in  r). 

The  change  in  +he  form  of  an  operator  is 


-  -  -  =  i[F,..], 

whe^e 

def 

F  -0G. 

Thus  the  time  development  of  a  system  may  be  visualized  as  a  sequence  of 
infinitesimal  unitary  transformations  with  the  generator  F.  F  is  the  time 
development  operator. 

If,  in  particular,  0^  =  C-  and  G  is  independent  of  time,  then  (taking  to=0), 


U  =  exp(  G-  ) 
and 

IT  ;  =  exp(  ..  -'G  ) *. 

Tne  operators  .  G  and  U  commute  and  therefore 


F  -  U  (  .  -  G  )  U”  *  =  w  G  . 

Then  defining 

def  (  •) 

A  (the  time  development  operator  if  the  symmetries  are 

.  -OM:n  '  one  can  write 

=i[  6-  ,  A  ]  =  i  .*••  (.'.* v)  f"  . 

t-u 


de  note  that 


Re- i (  • .  ) 


43 


which  implies 


ivic2(/\-v)‘f’-  =  -  ivK-'v2F-  . 


We  define 


def 

-•  ;  -  Tv'- 


then  , 


z  {.VK272f-2}  =  -  vn2.  (42 

t=0 

Equation  (42)  allows  a  simple  example  of  how  the  operator  algebra  can  be 
used  to  roughly  approximate  the  effects  of  the  interaction  on  the  structure 
decay  rate. 

One  assumes  that  ^M-i  remains  constant  for  times  of  the  order  of  an 

dt  IfO 

e-folding  time,  then  if  one  defines  t  as  the  e-folding  time  for  the  structure 
d(A  )  I 

dt  l.r  ’  vr>- 

■A  ( t)  • 


A?(t) |  the  change  in  A  (t)  in  time  t. 


however,  the  expression  for  ' 0 ( t ) >  is 


E ( t )  •  =  exp  •-  ■  (t-tQ)} |B(tQ)> 


if  the  symmetries  are  not  broken,  thus 


B ( t )  -  =  exp  (;■  +  h;)(t-t0)  | B ( tQ ) ' 


if  (t-t  )  =  ;  arid 


B( t)  =  exp  -  .(<•'  +  n2)(t-tQ)  j B( tQ) > 


if  the  symmetries  are  broken. 
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It  has  been  shown  elsewhere  (Norwood,  1969)  that  ,  where  L  is  of 

the  order  of  the  diameter  of  the  structure. 

One  can  now  make  a  rough  estimate  of  the  increase  in  the  decay  rate  of  the 
structures  interacting  at  the  center  cf  the  mirror  system  by  noting  that  f  is 
a  very  slowly  varying  function  of  (x,y,z),  thus 


f--  «  ■■■ 


This  implies  that  v~e'VK  ^ ,  i.e.,  the  decay  rate  is  not  altered  appreciably 
if  the  symmetry  is  slightly  broken  due  to  the  interaction. 

A  more  precise  estimate  can  be  easily  made  if  the  amount  of  symmetry  breaking 
can  be  estimated  more  accurately.  Returning  to  <  ~  f,  the  function  f(x,y,z)  can 
be  evaluated  if  the  expansions  for  *^(<51)  and  enNn($)  are  known.  If  these 
expansions  are  made  in  the  form  of  an  analytic  series,  the  amount  of  symmetry 
breaking  can  be  related  to  the  number  of  terms  retained  in  the  expansion.  For 
example 


def  ' 


with  a  Mn( : )  =  1  for  n  =  0 


*  M  ( ; ) ;  A  •  B  •  d  ■  |{1  -  V  P  ( (j )  }A*  Bd  r 

nn  !  n=0n 


and  P  U)  =  0  for  n  -  0, 


where  P  (:)  is  a  complete  set  of  functions  of  the  coordinates  whose  sum  has 
been  normalized  to  unity.  Then,  if  all  terms  from  n=0  to  n=<*  are  retained, 
the  ymmetry  correspondi ng  to  jA-Sdi  is  completely  broken.  If  only  the  n=0 
and  n- ,  terms  are  retained,  the  symmetry  is  weakly  broken,  etc.  Similar  ex¬ 
pansions  for  -■  N  ( ; ) B * vd :  and  l'*'  P  A*vdi  can  be  made.  A  more  complete  dis- 
nii  •  n  n 

cussion  of  these  expansions  will  be  given  in  another  paper. 


It  has  been  shown  elsewhere  that  each  of  the  integrals  can  be  associated 
with  one  of  the  independent  vector  fields  that  are  convected  with  the  flow 
(Wells  1970).  Thus  a  breakdown  in  the  convection  of  B,  v,  or  z  (the  generalized 
vorticity  vector)  can  be  estimated  from  experimental  observation  of  the  corres¬ 
ponding  interactions.  This,  in  turn,  will  enable  one  to  determine  the  number 
of  terms  to  be  retained  in  the  expansions  of  I-|  ,  and  1^,  i.e., 

't  -> 

;cl-|nMn(c){A-B}dT,  /BlnNnU)dv,  and  /YlnPnU)di. 

It  is  suggested  that  the  "equal  time  commutation  algebra"  represented  by 
Eq.  (42)  combined  with  a  knowledge  of  the  group  algebra  of  the  M  operators 

1.  V 

allows  calculation  of  some  detailed  properties  of  the  interactions  of  force- 
free,  coll  i  near  structures. 

The  difficulties  connected  with  operator  convergence  in  the  case  of  broken 
symmetries  for  quantized  operator  fields  (Orzalesi,  1970)  does  not  arise  in  this 
classical  theory.  Our  "current  algebra  approach"  is  only  an  analogy  employed  to 
clarify  the  classical  field  symmetries  involved  in  global  stability  calculations. 

V.  Conclusions 

Group  theoretical  and  symmetry  formalisms  have  teen  used  to  derive  and 
discuss  the  Chandrasekhar-Woltjer-Wentzel  variational  approach  to  finding 
equilibrium  solutions  of  the  conservation  equations  of  an  MHD  fluid  plasma 
model.  The  theory  has  been  derived  from  first  principles  utilizing  the 
multiple  integral  techniques  of  Caratheodory  and  Weyl .  The  principle  of  least 
constraint  is  shown  to  emerge  naturally  and  rigorously  from  these  formalisms. 

This  principle  combined  with  the  fundamental  variational  formula  for  multiple 
independent  integrals,  allows  calculation  of  the  global  stability  of  closed 
plasma  structures.  The  second  paper  of  this  set  presents  experimental  verifi¬ 
cation  of  the  most  important  predictions  of  this  theory. 
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Much  confusion  exists  in  the  current  literature  about  the  possible  values 
of  6  for  a  force-free  structure.  These  are  references  to  the  solution  of  the 
magnetohydrostatic  equation 

Vp  =  jxB.  (A  1) 

There  is  much  discussion  of  the  stability  of  plasma  configurations  which  are 
bounded  by  regions  of  vacuum  fields  for  which 

-*  ->■ 

jxB  =0  or  (vxB)xB  =  0. 

It  is  then  noted  that  one  must  have 


vp  =  0, 


and  since  this  condition  must  exist  in  a  boundary  where  pressure  gradients  must 
be  fini te,  then  p  0  and 


&  =  2=0- 


The  implication  is  that  force-free  fields  are  associated  with  low  6  plasmas.  Th 
is  true  if  tq.  (A  1 )  is  the  equilibrium  equation  for  the  region  of  the  plasma 

i 

under  discussion. 

The  plasma  we  are  discussing  in  this  paper  has  the  equilibrium  equation 

-*■  def. 


1 


0  -  -  v(p  +  w  ev2)  +  jxB  -  p(^xv) 


where 


vxv. 


1 


If  the  flow  is  force-free  and  coll  inear,  this  equation  takes  the  form 
•  (p  +  \  pv  )  =  0 

L. 

or  p  +  j  ,  v •  =  constant. 

If  the  flow  is  an  equi parti tion  flow  (the  most  stable  structure  has  g  =  — 
1  * 

and  v  =  — - —  B),  then  one  has 


1_ 

2 


,  V- 


and 

P  +  =  constant.  (A  2) 

o 

The  integral  of  over  the  surface  of  any  fluid  element  within  the  boundary 
Mo 

of  the  closed  plasma  structure  under  consideration  must  be  identically  zero 
since  the  body  force,  j>B,  is  zero  throughout  the  volume  of  the  structure. 

Tnus  Eq.  (A  2)  does  not  define  e  inside  the  boundaries  of  the  plasmoid.  The 
structure  must  be  supported  by  currents  somewhere  in  the  surrounding  plasma 
which  interact  with  the  vacuum  guidefield.  This  'leans  that  somewhere  on  the 
boundary  of  the  plasma,  presumably  close  to  the  walls  of  the  containing  vessel, 
there  must  be  a  transition  to  vacuum  magnetic  field  if  the  whole  configuration 
is  to  be  supported  free  of  the  walls.  Thus,  somewhere  near  the  walls,  e  must 
again  be  a  definable  quantit'  and  we  must  have  the  condition 


i 


This  means  that  the  configuration  consisting  of  the  plasmoid  and  its 
surrounding  plasma  is  limited  to  values  of  6  between  one  and  zero.  It  does  not 

mean  that  all  force- free  structures  in  a  dynamic  plasma  are  necessarily  low 
s  true  cures . 
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